Coordinate scaling in time-dependent current density functional theory 
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The coupling constant dependence is derived in time-dependent current density functional theory. 
The scaling relation can be used to check approximate functionals and in conjunction with the 
adiabatic connection formula to obtain the ground-state energy from the exchange-correlation kernel. 
The result for the uniform gas using the Vignale-Kohn approximation is deduced. 

PACS numbers: 31.15.Ew, 71.65. Gm 
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Time-dependent functional density functional theory 
(TDDFT) |l| is developing rapidly as a tool for predicting 
electronic response to laser fields, both weak and strong 
0. For weak fields, linear response applies, and perhaps 
the most popular present application of TDDFT is in 
calculating optical response of molecules, including tran- 
sition frequencies For strong fields, TDDFT allows 
prediction of many properties in response to intense laser 
pulses, such as high harmonic generation 0. It seems 
likely that TDDFT will play a key role in the emerging 
field of electron quantum control [3 • 

The time-dependent scheme can also be used for 
strictly ground-state properties by using the adiabatic 
connection formula to compute the static exchange- 
correlation energy H 1 Although using TDDFT 
to calculate ground-state properties might seem an 
unwarranted complication, the approximate exchange- 
correlation energy functional within this scheme has 
many useful properties. For example, it correctly de- 
scribes the static correlation for bond dissociation [8j and 
can be used to calculate accurately van der Waals disper- 
sion energies H G3 [H G3 . 

To use TDDFT in the adiabatic connection formula, 
one must generalize the response functions to arbitrary 
coupling constant A. In DFT, this has a very precise 
meaning, as the density is held fixed while the Coulomb 
repulsion between electrons is multiplied by A [E El- 
Coordinate-scaling is used to derive the A-dependence 
of quantities in DFT [T3|, or current DFT Q. It leads 
to many useful results, such as the virial theorem for the 
exchange-correlation energy 0] > an d exact conditions on 
that energy 0. It can be used to check that approxi- 
mate functionals have the right scaling behavior [lq. 

Timc-dcpcndent current density functional theory 
(TDCDFT) is a more general scheme where the cur- 
rent density is the basic variational parameter instead 
of the density and it can include arbitrary magnetic 
fields. Unlike TDDFT, TDCDFT can be approximated 
with local or semi-local functionals without any con- 
ceptual difficulty pj| and is now being used to calcu- 
late excitations of quantum wells 0, 0] , atoms [2(J , 
molecules Ell I22I I2H. 124) and single molecule trans- 
port |25l l2q . Five years ago, the connection between 
coordinate-scaling and the coupling constant was derived 
for TDDFT [27} • The present paper generates analogous 



results within the more general framework of TDCDFT. 

TDCDFT [HIH starts from the Schrodinger equation 
for N electrons in a vector potential ae X t: 
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(Pi + a cxt (f l , *)) + V cxt + V e( 



where V ex t 
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denotes the one-body potential, 
f j — r 7 - 1 1 and pi — — iVj. We use 
m e = h = 1), and there 



atomic units throughout (e 2 
is an implicit speed of light constant, c, included in 
the vector potential, i.e., B = cV x A, where B is 
the usual magnetic field. The physical results from 
the Schrodinger equation above are invariant under the 
gauge transformation 



v ext (r,t) = Vext(r,t) 



<9A(r,t) 
dt 



a ex t(r,0 = a oxt (r, t) + VA(r,t), 



(2) 
(3) 



where A is an arbitrary function. The gauge freedom can 
be used, for example, to remove of the scalar potential 
by choosing 9A(r, t)/dt = w cxt (r, t). 

As the density is the conjugate variable to the scalar 
potential v GXt (r,i), the conjugate variable to the vector 
potential is the current density 



1 N 

-^{^(^(r-rO + ^r-rOv^t)}, (4) 
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where v,; (t) = Pi + a ext (fi, t). The basic theorem of TD- 
CDFT [23, states that, for a given initial wavefunc- 
tion, a given j(r, t) is generated by at most one ae Xt (r, t), 
up to a gauge transformation. The density and the cur- 
rent are related through the continuity equation 



dn(r, t) 
dt 



•V-j(r,t) = 



(5) 



To derive the coupling constant dependence, we 
first transform the coordinates to (Ti,t) — (Ar^,A 2 £') in 
Eq. 0: 



iW^ + a oxt (ArUV) 



t 1 i d T 1 

W\ = W\. 

A 2 dt' 
(6) 
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where the prime means that the quantity is evaluated at 
(r',t'), V e ' x t,A = Z)i=i u ext(Ar^A 2 t'), and the scaled nor- 
malized wavefunction is <3/' A = X 3N / 2 *(Ar' 1 ...Ar' 7V , \ 2 t'). 
Consistent with Ref. j^, we define the scaled density 
by n\(r, t) = A 3 n(Ar, \ 2 t). Now we also define the scaled 
current density 



j A (r,i) = A 4 j(Ar,A 2 t). 



(7) 



Continuity (Eq. ©) remains satisfied for all A. Multi- 
plying Eq. © by A 2 and omitting the primes, 



j^S ( pl + >a C xt(Af l ,A 2 t)) 2 + A 2 K xt A + Al/ oc | * 



• 9 T 

(8) 



We define aZ ex t[j, ^o] as the vector potential for a sys- 
tem with modified coupling constant A, which gives rise 
to current j starting from wavefunction ^o- Thus we 
identify 

a ox t^[jA,%A](r,t) = Aa cxt [j,4<o](Ar,A 2 i). (9) 

Although the A-dependence of the external potentials 
in Eq. (JSJ) is generally complicated, by virtue of the 
one-to-one correspondence between current and poten- 
tials [2^, , the vector potential appearing in Eq. © 
is that unique potential producing current density j(r, t) 
from the initial wavefunction ^Sq, with electron-electron 
interaction AV^ C - 

Next we apply the same argument to the Kohn- 
Sham system, where the electrons are non-interacting 
(V cc — 0) and a oxt (r,i) is replaced by an effective vec- 
tor potential, a s (r,<), defined to reproduce the same 
current as the interacting system. Since our previous 
argument does not depend on the interaction, a^(r,t) 
also satisfies Eq. |J5}- And the Hartree vector poten- 
tial, a H (r,t) = V fdt' jdr' e 2 n(r',t')/\r - r'|, satisfies 
the same scaling. From the definition of the exchange- 



correlation potential, a s = a 



a xr , we see that 



it must obey the same scaling as the other vector poten- 
tials: 

a xc Z[j;*o,$o](r,t) = A a xc [j 1/A ; * ,i/A, $o,i/A](Ar, X 2 t), 

(10) 

where there is also a functional dependence on the initial 
Kohn-Sham wavefunction <&o, from a s (r,i). This is the 
central result of this work. 

When the vector potential is irrotational, i.e. can be 
gauge-transformed to a scalar potential, the TDDFT 
A-dependence of Ref. j2?| can be derived from these 
more general results. From the gauge transformation, 
Eqs. one can see that an irrotational vector poten- 

tial is transformed to a scalar potential through da/dt = 
Vu. Inserting Eq. (|1U[) . we find 

V*4(r,t) = A3a xc (Ar, \ 2 t)/dt = X 3 V x v xc (Xr, X 2 t), 

(11) 

where Va = d/d(Xr). Requiring the potential to vanish 
far from the system we recover w xc (r, t) — A 2 w xc (Ar, X 2 t) 



from Ref. This relation can also be derived directly 

from Eq. © by the same arguments used for the vector 
potential. 

While Eq. © represents the most general form, 
applicable to all TDCDFT applications, we next look 
at the special case of the linear response of an elec- 
tronic system. The susceptibility, x, is usually defined 
by Sn(r, t) — J dr'dt' x(r, t; r', t 1 ) Sv(r\ i'), where Sn is a 
small change in density due to a small perturbation in the 
potential, 5v. We sometimes represent the previous equa- 
tion as Sn — x*5v. Since we now have a vector potential, 
we can generalize the linear response to <5j = *x * We 
restrict ourselves to applying (time-dependent) pertur- 
bations on systems for which the external potentials are 
static. The response can then be considered as a func- 
tional of the ground-state density only, not the current. 

The scaling relation for the linear response exchange- 
correlation kernel in TDDFT is given in Ref. [3(J. In 

TDCDFT the tensor analog is defined as / xc = <5a xc /(5j 
and we can find the scaling relation with the functional 
differentiation 

a xc /[n + Sn](r,t) - a xc /[n](r, t) - 

A(a xc [n 1/A + <5n 1/A ](Ar, \ 2 t) - a xc [n 1/A ](Ar, \ 2 t)) = 

A J dr'dt' / xc [n 1/A ](Ar, r',A 2 t-t') <5ji/ A (r', t') = 

A 4 



A (X A dr)(\'dt) / xc [n 1/A ](Ar,Af, A (*-*)) 



Eq. (fTUfl implies 

/Qno](r,r',t-t') = A 2 ,fe[n , 1/A ] (Ar, Ar', \ 2 (t - i'ijl2) 

or, in frequency space, 

fhc[n ](r,r',w) = ,f xc [n 0il/A ](Ar, Ar',o;/A 2 ). (13) 

These results are needed to implement the TDCDFT ver- 
sion of the adiabatic connection formula as shown below. 
In the special case of a uniform electron gas 



/Zxc[no](q,w) 



- ho[n ,i/x] j$ 



(14) 



The above relation implies that, for a uniform gas, know- 
ing the exchange-correlation kernel as a functional of the 
density is the same as knowing the coupling constant de- 
pendence; this was used for the equivalent TDDFT case 

There have been various approximations proposed for 
/xc [H IH and IM E3 since they are such 
important quantities. The main TDCDFT approximate 
functional currently in use is the Vignale-Kohn (VK) 
functional |l7j |. This is the gradient expansion in the 
current density, and uses as input the q — > limit of 
both the longitudinal exchange-correlation kernel, /x C (w) 
(which is precisely the scalar / xc (w) of TDDFT), and the 
transverse kernel, / t c (cj) of the uniform gas. We have 
checked that the VK functional respects the above scal- 
ing relation, Eq. I|13f) . provided that /ic' T ^( w ) used in 
constructing the functional also respect the appropriate 
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scaling. The most recent approximation for these ker- 
nel components is that of Qian and Vignale [37]. We 
verified that it satisfies Eq. 1)14(1 . assuming the Landau 
parameters are invariant under simultaneous scaling of 
the density and the coupling constant. 

Just as for the exchange-correlation potential, the scal- 
ing relation for the exchange-correlation can be derived 
from TDCDFT. When the vector potential is irrota- 
tional, the scalin g re lation of / xc reduces to that of the 
scalar kernel / xc [23, via 

VV7xoH(r,r» = W 2 /Qn ](r,r', W ) 

= w 2 /xcKi/A](Ar,Ar>/A 2 ) 

= A 4 V A v^/ xc [ no i/A ](Ar,Ar',^/(^) 

Then, since / xc -> as r -» oo for any finite sys- 
tem, integration implies than the scaling of the cur- 
rent kernel reduces to the scaling of the scalar ker- 
nel, /4jn ](r,r',w) = A 2 /xc[n ,i/A](Ar, Ar', uj/X 2 ), as in 
Ref. 30]. 

Similarly to the adiabatic connection formula used 
in ground-state DFT H E H|], which relates the 
exchange-correlation energy to the susceptibility, we in- 
troduce the adiabatic connection for the ground-state of 
a system with a static scalar potential using current DFT 
susceptibility 

E XC = -i d X £ |i Tr (P - M r) l] * (16) 

the trace is Tr (a) — J dr^2 i a„(r, r) and 

= -W cc V/w 2 . The symbol 1 stands for S(r-r') 
The tensor susceptibility is related to the exchange- 
correlation kernel through 39] 

Xl = X? + X? (AT + fVc) 7l , (17) 

where is the tensor susceptibility for the Kohn-Sham 
system 

X s,ij(r,r',uj) = n (r) 1 + 

f , ^(r)V 8 ^(r)^(r')V^ a (y;^ 

where / is the occupation number, i.e. 1 for an occupied 
state, for an unoccupied one, and r\ is infinitesimal. 
The Kohn-Sham wavefunctions and energies are denoted 
by (f> and e. 

In the special case of a homogeneous gas, the longitu- 
dinal and transverse responses decouple, and reordering 



the terms within the trace of Eq. (|16|l shows that only 
the longitudinal components contribute to £y C , i.e., it 
reduces to the usual scalar case. Lein et al |23 tested a 
variety of approximations to the scalar / xc for the uni- 
form gas, to see how well they reproduced the known 
correlation energy. To perform the same test for the VK 
functional, we first note that, although VK is a gradient 
expansion in the current, yielding terms of order q 2 , these 
terms are actually zero-order in q when transformed back 
to the equivalent scalar kernel via Eq. 1)15(1 . So we find 
that VK, inserted in the current adiabatic connection 
formula, reduces to inserting f xc {uj) — / x £ lf ((Z — > 0,w) 
in the usual scalar adiabatic connection formula. This 
approximation was already tested by Lein et al, and is 
labelled 'local RA' in their work. (Although they used a 
different parametrization [35j from QV |3jfl , the results 
are unlikely to depend strongly on such details.) They 
found about a factor of 2 reduction in error relative to the 
adiabatic local density approximation (ALDA). We have 
thus demonstrated that, for the special case of the uni- 
form gas, the VK approximation, inserted in the current 
adiabatic connection formula, improves over ALDA. 

Carrying out a calculation of Eq. i(16|) on molecules 
or solids is much more computationally demanding than 
the usual ground-state calculations with approximate 
exchange-correlation energy functionals, but is probably 
not much more expensive than the scalar case. Such cal- 
culations are presently being performed [H, |4jj because 
the use of the adiabatic connection formula correctly de- 
scribes the dissociation of molecules [8j and dispersion 
energies 0, E3 ■ The exchange-correlation kernel 

of TDCDFT being better suited to local or semi-local ap- 
proximations than the pure density theory |l7j ] , we would 
expect that it would supersede TDDFT when used within 
the adiabatic connection formula. 

To summarize, we have used coordinate scaling to de- 
rive the coupling-constant dependence of the exchange- 
correlation potential in TDCDFT. We have derived the 
adiabatic connection formula for TDCDFT, and shown 
how the VK approximation performs for a uniform gas. 
We have also given explicit formulas relating both po- 
tentials and kernels in TDCDFT to their couterparts in 
TDDFT. Given both the recent use of TDDFT exchange- 
correlation kernels in the adiabatic connection formula, 
for calculating bond dissociation curves, and the use and 
tests of TDCDFT for excitations in which TDDFT has 
shown limitations, it is clear that an important appli- 
cation of this work is likely to be realized in the near 
future. 
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